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Abstract 



We critically examine the evaluation of the hadronic contribution to the running of 
ON ' 

the QED coupling, a, from Q 2 = to Q 2 = M|. Using data for e + e~ — ► hadrons we 



find that ce(Mf ) _i = 128.99 ±0.06, as compared to the existing value of 128.87 ±0.12. 



*On leave from the Department of Physics, University of Durham, England 



The improvement in the measurements of electroweak quantities allows high precision 
tests of the Standard Model in which the measured Z boson mass is related to other ob- 
servables (see, for example, the recent reviews in Refs. |l], El). Surprisingly, out of the 
three accurately measured quantities (a, Gf, Mz) which determine the Standard Model, 
the largest uncertainty comes from the running of a from Q 2 = 0, where it is precisely 
known, up to the Z pole, which is the scale relevant for the electroweak precision tests. In- 
deed, other electroweak quantities are being measured with an accuracy comparable to that 
associated with a(M z ). The source of the ambiguity in the value of a(M z ) is the hadronic 
contribution to the photon vacuum polarization H(s). This contribution is determined by 
the dispersive sum of all possible hadronic states produced in e + e~ annihilation into hadrons 
via an intermediate photon 

as„ r°° R^(s') 



with a 1 = 137.036 and 



RelL^) = —P v ' ds' } (1) 

3n J Ami s \ s — s ) 

R11 = a(e+e~ -> 7 -> hadrons) 

Here a^ = 47ra(s) 2 /3s is the lowest order point-like e + e~ — > 7 — > fi + fi~ cross section. The 
fifi cross section is expressed in terms of the running coupling a(s) in order to eliminate any 
QED effects from the hadronic contribution to the current-current two-point function Hh/oc. 
The currently accepted determination of the running of a is based on the analysis of 
Burkhardt et al. jl. They used data for i? 77 , supplemented by narrow resonance contribu- 
tions, to obtain a hadronic contribution to the running of a of 

Aa h (M 2 z ) = -Re U h (M z ) = 0.0288 ± 0.0009 . (3) 

When the leptonic contribution, 



rv 

Aa,(Mi) = f £ 
on -. 



m^M + o( 



m. 



mj 3 v Mf' 



0.03142 , (4) 



is added, their result gives 

- Ren(Mf) = Aa = Aa fc (M|) + Aa e (M z ) = 0.0602 ± 0.0009 , (5) 

which translates to 

a(Ml)- 1 = (1 - Aaja" 1 = 128.78 ± 0.12 . (6) 



The analysis 0] was subsequently updated by Jegerlehner |jj to give 

Aa h (Ml) = 0.0282 ±0.0009, (7) 

a(Mf)- 1 = 128.87 ±0.12. (8) 

Here we are using the effective QED coupling, which is denoted by a in the Review of Particle 
Properties [|T]. Q. 

For electroweak precision tests it is important to see if the determination of ReII/ l (M|) 
(and hence of a (Mf )) can be improved and, in particular, the error reduced. In the following 
we compare our analysis with the original work of Ref. Q since it lists the various contri- 
butions to Aah in detail. Indeed the first column of Table fjj which is taken directly from 
Ref. [|J , shows the contributions to Re 11^ from different W = v? regions of the dispersion 
integral, together with the associated errors. We note that the largest error arises from the 
region 2.3 < W < 9 GeV, which contains the cc resonance region together with the 'contin- 
uum' regions below both the cc and the bb threshold. The evaluation of the contributions 
from this region, used in Ref. || f\, relied on the original MARK I |7| data for R(W). In fact 
the major uncertainty in the determination of Re IT/j is associated with the normalization 
errors of the measurements of er(e + e~ — > hadrons). Thus the 10% error associated with the 
2.3 < W < 9 GeV contribution of Ref. M reflects the 10% normalization uncertainty of the 
MARK I data. 

In the continuum regions well above the qq thresholds we are now in a position to use 
perturbative QCD to predict R 11 extremely accurately. First our knowledge of the QCD 
coupling has considerably improved since the previous calculation Q of Reil^. Even if 
we take a conservative view of the average of all of the measurements of a s M, we can 
conclude a s (M|) = 0.118 ± 0.007. Moreover i? 77 is known up to, and including, the 0(al) 
contributions, and the running of a s is known to 3 loops. Finally we now have experimental 
evidence of the value of the top quark mass m t . Thus, for example, if we evaluate R 11 at 
values of W = 3, 9 and 150 GeV just below the cc, bb and ti thresholds respectively, we find 



tThe updated analysis pi uses data from the Crystal Ball collaboration 



i? 77 = 2.17 ±0.03 atW = 3GeV, (9) 

= 3.58 ± 0.04 at W = 9 GeV , (10) 

= 3.80 ±0.01 at IV = 150 GeV. (11) 

We allow for the change in the number of flavors at each qq threshold both in a s (s) and 
in i? 77 (s). The errors include the ±0.007 uncertainty in the input value of a s (M§), an 
uncertainty from yet unknown 0{oQ) terms (taken to be equal in size to the a^ contribution 
to i? 77 ), and uncertainties from threshold effects. The main threshold uncertainties arise 
in the bb channel because we combine data on the resonance contributions with the QCD 
formula. We estimate these effects by comparing a naive (3(3 — /3 2 )/2 threshold behaviour 
with the full 0(a s ) QCD formula || and find that bb threshold uncertainties contribute very 
little to the error on a(M§). 

Given that i? 77 is known so precisely in the continuum regions, we may use it to improve 
the normalization of the experimental measurements of a(e + e~ — > hadrons). Such a program 
was carried out for 21 experiments by Marshall |J in a detailed study performed in 1988. 
He noted that many experiments which partially overlap the MARK I region have smaller 
systematic errors than the MARK I data. As a result of a global QCD fit to the world data 
for R he concluded that the experimental normalization of 1.00 ± 0.10 of the MARK I data 
should be adjusted to 0.850 ± 0.019 to bring them into line with the world data (which 
was by far the biggest adjustment of data that he obtained). In their paper [|7| MARK I 
quote a normalization error of ±20% at W = 2.6 GeV decreasing smoothly to ±10% for 
W > 6 GeV. Marshall's adjustment brought the MARK I data into excellent agreement with 
QCD expectations in the continuum regions below the bb and cc thresholds. Clearly such an 
adjustment will have a dramatic effect on the value, and the accuracy, of a (Mf ). Incidentally 
at the time of Marshall's analysis the coefficient of the (a s /7i) 3 term in i? 77 was erroneously 
large and of the wrong sign (±65 instead of —12.8 for five flavors). As a consequence the 
renormalization factors found by Marshall should have been even smaller. 

To evaluate (|l|) we assume that R" 11 is given by perturbative QCD in the continuum re- 
gions 3 < W < 3.9 GeV and 6.5 < W < oo, apart from the if; and T resonance contributions. 



Typical errors on R 11 in these regions are shown in Eqs. ([|)-(11). Following Marshall's pro- 
cedure, we also use the continuum values of i? 77 to normalize the various data sets. For the 
MARK I data [[/]] we find that an overall renormalization of 0.83 ± 0.02 is required. In fact 
fitting to the MARK I data in the W < 3.85 GeV continuum region gives essentially the same 
renormalization factor as fitting to their W > 6.5 GeV continuum data, see Fig. |T](a). One 
option would be to use the renormalized MARK I data to evaluate the contribution to (|l|) 
from the interval 3.9 < W < 6.5 GeV, between the two continuum regions. To be precise we 
could use the curve of Fig. |l|(a) in which the portion between 4.5 < W < 6.5 GeV is shown 
dotted. To check this 'MARK I' curve, we compare with the more recent and precise Crystal 



Ball measurements |6|, 10|. First we renormalize the Crystal Ball ('90) measurements us- 
ing their W > 6.5 GeV data. A factor 1.06 ±0.02 is found. Indeed with this renormalization 
all the Crystal Ball ('90) data are well described by i?(QCD), see Figs. |](a,b). From Fig. |l|(a) 
we see that the more precise Crystal Ball ('90) data lie significantly above the renormalized 
MARK I data for W ~ 5 GeV. Evidence in favor of the higher Crystal Ball values and, in 
particular, for taking the continuous (rather than the dotted) curve for R, comes from two 
first-generation experiments, DASP (Tl| and PLUTO fL2f , see Figs. |](c,d) respectively. The 
data are shown after they have been renormalized [J] to i?(QCD) at W ~ 3.6 GeV. On the 
other hand for W < 4.5 GeV the data of Figs. [l|(b,c,d) show that the line drawn through 
the MARK I data is a reasonable representation of R in the cc resonance region. We note 
that both the Crystal Ball ('86) |Uj and the DASP (JTTJ data are able to resolve the ^(4.04) 



and ^(4.16) states and, especially, that our curve is a good average of R for this resonance 
region. The contributions to Kellh, obtained from integrating R over the continuous curve 
in Fig. [I], are shown in Tables [I] and ^|, together with the contributions of the families of ijj 
and T resonances. The resonance contributions are determined from [|J] 

n res = -^ a 2 , (12) 

M a (M 2 ) 2 ' l ; 

where M and T ee are the mass and leptonic width of the resonance, respectively, and where we 
account for the running of the effective QED coupling at the resonance scale. Equation ( |12"D 

follows from integration over a narrow Breit-Wigner resonance form. 

*We obtain similar normalization factors to Marshall H, after allowing for the change in the 0{a 3 s ) 
contribution. 



The errors associated with the continuum contributions are estimated as for Eqs. (||)- 
(|TT|). For the intervening interval, 3.9 < W < 6.5 GeV, we estimate the error by repeating 
the calculation using the dotted line in Fig. |l](a). The contribution reduces from 2.90 x 10~ 3 
to 2.71 x 10~ 3 , and we regard this change as representative of the uncertainty of this interval. 

We now turn to the region below W = 3 GeV. Here there are many experiments measuring 
e + e~ annihilation to specific hadronic channels. We evaluate the contribution to Re 11^ from 
this region in four separate parts; see Table 0. First we calculate the contribution from 
e + e~ — * 7r + 7r~ by integrating over R obtained from detailed measurements []n| of the pion 
form factor F n (s) via 

fl(a) = 4 (l - -y*J \F n (s)\ 2 ■ (13) 

To be specific we integrate over the curve shown in Fig. || and assign to this contribution 
a ±4% uncertainty arising primarily from the experimental normalization of l-F^s)! 2 . This 
region is dominated by the p resonance (with the resonance shape mutilated by p-u interfer- 
ence). For example the intervals 1 < W < 1.4 and 1.4 < W < 2 GeV only give contributions 
to —Re lift of 0.16 x 10~ 3 and 0.02 x 10 -3 respectively. Secondly, we include the contribu- 
tion due to the cj(782) resonance using (|T2j), where the error reflects the uncertainty in the 
observed leptonic width ||. Thirdly, we calculate the e + e~ — » KK contribution to RelT/j 



using the parametric form for the kaon form factor determined by Bisello et al. |14| from a 



fit to e + e — ► K + K data |T4L |T5|. By far the dominant contribution here comes from the 



<p resonance, though there are small contributions from the p, u — > KK resonance tails and 
an even smaller contribution from the 1.5-1.6 GeV resonance region; see the dotted curve 
in Fig. |3|. The error reflects the uncertainty in the <fi leptonic width [Q. To a good approx- 
imation the total contribution is twice the K + K~ contribution. On the resonance this 
allows for K°K° and 3ir contributions, while above the resonance it represents the total KK 
contribution since, away from threshold, the effect of the K + K~ — K°K° mass difference is 
suppressed. Finally we have the contributions to Rell/j from the region above W = 0.9 GeV 
due to multi (> 3) pion production. For the region 0.9 < W < 1.45 GeV we use the sum of 
the data for the exclusive channels ||16|| . The dominant contribution comes from 7r + 7r~7r 7r° 
and 7r + 7r~7r + 7r~ production. Above W = 1.45 GeV we use a line through the R data of 
Refs. |nj, joining on to the QCD value at W = 3 GeV, as shown in Fig. [| 



It is instructive to compare our results with those of Ref. 0. From Table [I] we see that 
the main improvement in accuracy comes from the 2.3 < W < 9 GeV and 12 < W < oo 
regions, due mainly to our use of i? 77 (QCD). The difference in size of the 2.3 < W < 9 GeV 
contribution can be attributed to our use of renormalized MARK I and Crystal Ball '90 
data. To make an exact comparison with Ref. [Q we should let m t —>■ oo, rather than taking 
the value m t = 174 GeV that we have used to include the e + e~ — > tt contribution. The 
resulting effect is that — ReII/ l (M|) would increase very slightly to 27.39 x 10 -3 . 

Although we have numerically integrated over p and <ft resonant shapes, and carefully 
considered individual contributions to e + e~ production processes, we find that there is a 
comparatively modest reduction in the error of the contribution to RelT/j from the low 
energy, W < 3 GeV, region. From the first five rows of Table [TJ we find that the contribution 
for W < 2.3 GeV is -(6.06 ± 0.25) x 10~ 3 as compared to -(6.34 ± 0.43) x 10~ 3 of the 
previous calculation 0. In fact we see from Table |2] that the error ±0.33 x 10~ 3 on the 
W < 3 GeV contribution limits the present accuracy of a(M§ ). 

In conclusion we find that 

a(Mf)- 1 = 128.99 ±0.06. (14) 

The large reduction in the error, in comparison to that found in Refs. Q, [| , arises because we 
use the precise perturbative QCD prediction for K 11 in the regions well above qq thresholds. 
As a consequence we find (i) that the contribution from the high W = y^ region is well 
determined, and (ii) that we are able to reliably normalize the data in the cc resonance region 
which significantly reduces the value and the uncertainty of a(Mf ). Although perturbative 
QCD stabilizes the contribution from the W > 3 GeV region, the contribution from lower 
values of W relies directly on the available data for e + e~ annihilation into hadrons. In 
fact the dominant uncertainty arises from the 1 < W < 3 GeV region and, in view of the 
importance of an accurate value of a(M|) for precision tests of the Standard Model, it is 
crucial to improve the accuracy of the low energy measurements of a(e + e~ — > hadrons). 
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Table 1: Contributions 


to -1000ReiI h (Mf). 




W range (GeV) 


Burkhardt et al. |§ 


This work 


2m^-1.0 




— 




3.21 


± 


0.14 


P 


3.484 


± 


0.171 




— 




us 


0.347 


± 


0.021 


0.31 


± 


0.01 


1.0-2.3 


1.981 


± 


0.391 


1.95 


± 


0.21* 





0.528 


± 


0.024 


0.59 


± 


0.02* 


2.3-9.0 


7.208 


± 


0.721 


6.56 


± 


0.22* 


V>'s 


1.084 


± 


0.057 


0.92 


± 


0.05 


9.0-12.0 


1.686 


± 


0.169 


1.67 


± 


0.06* 


T's 


0.118 


± 


0.005 


0.12 


± 


0.01 


12.0-oo 


12.368 


± 


0.371 


11.99 


± 


0.09* 


Total 


28.8 


± 


0.9 tt 


27.32 


± 


0.42 


a(Ml)- 1 


128.78 


± 


0.12** 


128.99 


± 


0.06 



* Includes the full e + e~ — > KK contribution. 

*These errors have to be added linearly. 

*Part of these errors have to added linearly and a part in quadrature. 

t1 ~An updated analysis § gives -1000Ren h (M|) = 28.2 ± 0.9 and 

a(Mf)- 1 = 128.87 ±0.12. 



Table 2: 


Summary of the evaluation of ReIl/ l (M|) 






W range (GeV) 


Information used 


-1000 Re n^Mf) 


2m n -3.0 










(i) 7T + 7T- (W < 2) 


Data for \F n \ 2 


3.39 


± 


0.14 


(ii) uj — ► 37r 


Narrow resonance formula 


0.31 


± 


0.01 


(iii) if.K" inc. 


Parametrization of e + e _ — »■ if K data 


0.59 


± 


0.02t 


(iv) e + e~ — > hadrons 
excluding 7T7T, if-K" 


i?(data), e + e~ — ► 47r etc. 


2.72 


± 


0.30 


3.0-3.9 


^(QCD) 


0.88 


± 


0.01 




+JM V', ^(3.77) 


0.92 


± 


0.05 


3.9-6.5 


Renormalized i?(data) 


2.90 


± 


0.19 


6.5-oo 


R(QCD) 


15.49 


± 


0.15 




+T(nS), n = 1, . . . , 6 


0.12 


± 


0.01 


Total 




27.32 


± 


0.42 



^Includes the — ► 3tt contribution. 

The errors are added in quadrature, except those for i?(QCD). 
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Figure Captions 



Figure 1: The (a) MARK I [§, (a,b) Crystal Ball '90 @, (c) DASP |T|, and (d) PLUTO p 
data for i?, respectively renormalized by 0.83, 1.06, 0.94, and 0.96, - factors which are found 
by fitting to i? 77 (QCD) in the continuum regions W > 6.5 GeV and W < 3.9 GeV. The 
Crystal Ball '86 data |10| in (b) are not renormalized. The continuous curve is the same in 



each plot and is the representation of the data used to calculate a(Mf ). The dotted curve 
is a representation of the MARK I data for 4.5 < W < 6.5 GeV which is used to estimate 
the uncertainty. The J/ip, if)' and if)(3.77) contributions are included using (|12|). 



Figure 2: Data on the pion form factor [13 . The curve is the representation of the data 



which is used to evaluate the rnr contribution. 

Figure 3: The dashed and dotted curves are a representation of the data for the contribution 
to R from e + e~ — > 7r + 7r~ and KK respectively. The continuous curve represents the data 
for multi (> 3) pion production. For W < 1.45 GeV it represents the sum of the multipion 
exclusive channels. Above W = 1.45 GeV the multipion data are taken from Refs. 0, 1? 



Also shown are three MARK I measurements 0, renormalized by a factor of 0.83. We 
numerically integrate over the p and 4> resonance forms, but include the a; (782) contribution 
via fin. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411377v3 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411377v3 
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